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Abstract 

An abstract formulation of quantum dynamics in the presence of a 
general set of quantum constraints is developed. Our constructive pro- 
cedure is such that the relevant projection operator onto the physical 
Hilbert space is obtained with a single, common integration procedure 
over the original Lagrange multiplier variables that is completely in- 
dependent of the general nature of the constraints. In the associated 
lattice-limit formulation it is demonstrated that expansion of the con- 
straint operator contribution to second order in the lattice spacing is 
necessary while, as usual, only a first-order expansion is needed for the 
dynamical operator contribution. Among various possibilities, coher- 
ent state path integrals are used to illustrate a completely functional 
representation of the abstract quantization procedure. 

1 Introduction 

The abstract operator formulation of quantum mechanics is generally recog- 
nized as being the most fundamental, while those formulations that involve 
one or another representation need to contend with representation-dependent 
issues (e.g., lack of global coordinates) that have no place in the essential 
formulation. Our goal here is to study quantum dynamics and quantum con- 
straints as much as possible in an abstract, representation-independent way. 



1 



Only then should one feel comfortable about introducing a representation for 
further analysis. Let us start with some basics. 

The evolution operator, from time t — to time t = T > 0, for a general 
quantum system with an abstract, self-adjoint, time-independent Hamilto- 
nian operator Ti, is given, in units where H — 1 (which are often used), by the 
unitary operator e~ t7iT acting on an abstract Hilbert space known hence- 
forth as the "original" Hilbert space. If the Hamiltonian is time dependent 
and self adjoint for alH, < t < T, the evolution operator is instead given 
by the unitary operator 

where T denotes the time-ordering operator. We may append to H(t) a sum 
of additional, time-independent, self-adjoint operators {& a }a=i with general, 
real- valued, time-dependent coefficients {\ a (t)}^ =1 leading, by assumption, 
for all t, < t < T, to the self-adjoint operator T~C(t) + \ a (t)Q a , with sum- 
mation implied (the generalization to a set of time-dependent, self-adjoint 
operators {® a (t)}a=i i s briefly discussed later). Taken as the generator of 
time translations, the extended evolution operator reads 

Te -iJ'llH(t)+\ a m a ]dt _ 

If the functions {A a (t)}^ =1 are fixed from the beginning and held fixed, the 
term \ a (t)§ a provides an addition to the Hamiltonian and modifies the dy- 
namics in the expected way. If we regard the functions {\ a (t)}^ =1 as variable 
external sources, then functional derivatives with respect to them lead to pro- 
totype expressions from which time-ordered correlation functions may be ob- 
tained. Alternatively, if the functions {A a (t)}^ =1 denote suitable stochastic 
variables, then an expectation of the given expression may relate to systems 
interacting with random potentials. Besides these examples, let us consider 
a rather different interpretation and application of the extra term \ a (t)Q a , 
which we now describe. 

Let us accept the hypothesis that the set {$ q }q =1 consists of constraint 
operators, which, when appropriately satisfied, direct attention to a subspace 
of ft called the "physical" Hilbert space and which is denoted by i3 p h ys - As a 
closed subspace, it follows that i3 p h ys = for some unique projection oper- 
ator E (= = E 2 ). In turn, the projection operator E should be given by 
an acceptable function of the constraint operators {$ q }q =1 , and it is the role 
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of an integration over the variables {\ a (t)}^ =1 to build the necessary projec- 
tion operator E. Hereafter, we shall refer to the set of functions {\ a (t)}^ =l 
as the Lagrange multiplier variables. 

In simple cases it is easy to determine a suitable projection operator E 
Hfl. For example, if A = 1, and so there exists a single constraint, then — see 
Appendix A — we observe, for any 5 = 5(h) > (N.B. 5(h) is not a Dirac 
delta function but a rather general function of h), that 



If $ = is part of a purely discrete spectrum near zero, then it is possible 
to choose 5 > such that E ($ 2 < 5 2 ) = E ($ = 0). When $ = is part 
of the continuum, we choose < 5 <C 1; one procedure to take the limit 
5 — > to arrive at the true physical Hilbert space ijtmephys is discussed 
in Appendix B. From this discussion, it follows in the continuum case that 
we may regard i3 p hy S as a provisional physical Hilbert space, which, despite 
the qualifying adjective, is, for an extremely small 5, already acceptable for 
calculational purposes. As a second example consider that the constraints are 
the generators of a compact Lie group, and let dm(\) denote the normalized 
group invariant measure, J dm(X) — 1. In that case (e.g., |], B HI) 



In order to generate the projection operator E , these examples illustrate 
that different measures for the Lagrange multiplier variables may be involved. 
Indeed, in order to construct E , even for a single set of constraints, it is 
generally the case that the measure for the Lagrange multiplier variables is 
not unique. 

The use of different measures for different situations is not wrong but 
it does mean that some specific properties about the system, especially the 
constraints, need to be determined before an expression for E can be found. 
This is clearly a nuisance, and it raises the question whether there exists a 
single, common, integration procedure by which we can construct E , and 
which, moreover, is a procedure that is completely independent of the specific 
constraints at hand. The answer to this question is "Yes!", and the aim of 
the present paper is to develop a universal prescription for determining E 





dm(\) = E ($ Q = 0, 1 < a < A) = E ($ a # 



^ = 0) . 
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from a given set of constraint operators, whatever kind of constraints they 
may be, by means of a prescription that is also independent of the specific 
Hamiltonian operator. 

There are various kinds of constraints that need to be considered. In the 
usual nomenclature (see, e.g., Q), the constraints may be: first class, closed 
or open; anomalous or second class; primary or secondary; irreducible or 
reducible; regular or irregular. In the reducible case, where the constraints 
are linearly dependent, it is possible that there is an infinite degeneracy when 
A = oo even if the number of degrees of freedom is finite; see below for an 
example. Regularity refers to whether the classical constraint is taken (say) 
as q = or as q 3 = 0, etc., along with the quantum repercussions of this 
kind of multiple choice test. Implicitly, we assume that the Hamiltonian and 
the constraints satisfy some minimal consistency condition, in particular, 
that the set of constraints already includes all the necessary constraints for 
a given Hamiltonian.. 

In what follows, all varieties of constraints will be dealt with in a uniform 
fashion. 



2 Construction of the Physical Space 
Projection and Evolution Operators 

We first observe that the extended evolution operator may be written in the 
form of a lattice limit given by 



lim Yl [Te '"' 1 1 

Kn<N 



where e = T/N and the directed product (symbol < — ) also respects the 
time ordering. Thus, this expression is simply an alternating sequence of 
short-time evolutions, first by A a (t)$ a , second by TC(t), a pattern which 
is then repeated N — 1 more times. The validity of this Trotter-product 
form follows whenever 7~C{t) 2 + $ Q 5 a/3 $# is essentially self adjoint for all t, 
< t < T. As a slight generalization, we shall assume that 7i(t) 2 + < l ) Q M a/3 $ / g 
is essentially self adjoint for all t, < t < T. Here the real, symmetric 
coefficients M Q/3 (= M Pa ) are the elements of a positive-definite matrix, i.e., 
{M a/3 } > 0. For a finite number of constraints, A < oo, it is sufficient 
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to assume that M a/3 = 5 alS ; other choices for M alS may be relevant when 
A = oo. 

With all this in mind, we shall explain the construction of a formal inte- 
gration procedure whereby 

J Te~^ (-^ Xamadt VR(X) =E($ a M Q % < 5(h) 2 ) , 

and for which the integral represented by J ■ ■ ■ T>R(X) is independent of the 
set of operators {$ Q } and the Hamiltonian operator TL(t) for all t. We will 
define the left-hand side of this expression by means of an additional lattice 
limit. 



Constructing the basic projection operator 

For any given n, 1 < n < N, let us divide the time interval ((n — l)e, ne] into 
N', N' < oo, steps of equal duration, e' = e/N' = T/(NN') > 0. In this way 
we are led to consider the integrated directed product 



J] e~ ie A "(«.o* a dR(X r (n^) 

1<1<N' 

for some measure dR(X r ( n j)) to be determined. Observe that in the indicated 
expression the continuum variable t has been replaced by steps of size e' 
consecutively numbered by r(n, I) = (n — 1)N' + (/ — 1) for 1 < n < N and 
1 < I < N'. As far as the previous integral goes, the particular value of n 
is unimportant and we will suppress it for the rest of this discussion. Thus 
for the present we replace r(n, I) by / for simplicity. Let us proceed in small 
steps. 

Introduce a measure dS{X{) for each A; and an integral over A = {A/}^ 1 
given by 

J= f J] e^'^dSiXt) 

J 1<1<N' 

= / II [l-^A^ a -ie /2 ArAf$ Q ^ + O(e /3 )]^(A0. 

J 1<KN> 
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We choose dS(Xi) as a normalized probability measure with the first two cor- 
relation functions (moments) given by JXfdS(Xi) = and JXfXf dS(Xi) = 
(27/V) M Q/3 , 7 > 0, values that are independent of the time-slice label /. 
With a minimal requirement on higher moments (such as finite fourth-order 
moments), it follows that 

j= n [i-e'7(^^%)i+o(^) 

1<1<N' 

where the last expression follows because the factors in the directed product 
are now identical for every time slice I. At this point we may take the limit 
N' — > 00, e' — > with N'e' = e held fixed. Stated formally, the net result is 
that 

a - vfo!-i)« Aa w* a *T><7m - P -«7(*«^*ri 



independently of n, 1 < ra < N. Note that many choices of the basic measure 
dS(Xi) will lead to the same result. The mechanism which lies behind this 
conclusion is similar to one leading to the central limit theorem in probability. 

The next step in the process is to analytically extend the variable 7, 
7 — > 7 + ir, where 7 > still and —00 < r < 00, namely 

e - £1 (*„M [ " i ^) _^ e - e ( 7 +ir)($ a M«' 3 $ / 3) 

and after achieving this extension, the following step is to pass to the limit 
7 — > + . Both of these steps are achieved by the single integral operation 

J e - E #„M^ rf(|r(7) 

= lim ^Bj2^e\ BT ' 2 / ^r-^T 2 e -e7(^M<-^ fl ) ^ 

The final step in the construction — see Appendix A — involves an integration 
over r given by 



e --(W^) rfT (r) 
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lim lim f L e-^M^M^ + Qr}^ 

;^0+ L->cc J_ L TTT 



C 

= E(e$ Q M a/3 $ /3 < e5 2 ) 
= E ($ a M a/3 $/3 < <5 2 ) , 

which achieves our goal. We occasionally symbolize the collective set of 
operations by J ■ ■ • VR(\), leaving the integrals over 7 and r implicit. 

To summarize, let us assemble the several, separate operations together. 
In particular, we observe that 



E($ a M a/3 $ /3 < 5 2 ) = J VR{\) Te 



. \ a (t)$ a dt 

(n— l)e 



sKm , im f L M^ + Qr] dT 

^-+0+ L~*oo J_ L TTT 

coo 

X 



B 

X 



poo 

lim ^/B/2^e^ Br2 / e^"^ 2 d 7 
I VS(\)Te<^ Xam « dt , 



where the proper definition of J ■ ■ ■ T>S(X) has been given above. In accord 
with the discussion in Appendix A, we note that if the final limit is replaced 
by lim f _>()-, the result becomes E($ Q ,M a/3 $ /3 < 5 2 ). 

A few comments are in order. It is noteworthy that J VR(\) = 1. Thus, 
in the hypothetical case that all the constraint operators are bogus, i.e., 
$q, = for all a (i.e., not simply weakly zero but strongly zero), it would 
follow that 

J 1VR(X) =E(0 < 5 2 ) = 1 . 

If some, but not all, of the constraint operators are bogus, then, even though 
an integration over the associated Lagrange multiplier variables takes place, 
the answer will be identical to one in which the Lagrange multiplier variables 
for the bogus constraint operators were ignored from the very beginning. 
As another hypothetical case, assume that there is an infinite number of 

identical constraint operators, for example, $ a = $ for all a, a — 1, 2, 3, 

Consequently, the extended evolution operator involves ri(i) + \ a (t)§ a = 
ri(t) + S^ =1 A a (t) $. We clearly see that there is only one relevant Lagrange 
multiplier variable, but to act on that knowledge would amount to choosing 
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an irreducible set of constraints from the original set. Instead, let us keep 
all the Lagrange multiplier variables as independent quantities and proceed 
along the lines developed in this paper. In this case let us choose, for example, 
M a>3 = {Q/n 2 a 2 ) 5 a/3 . Then it follows that the result of the integration with 
respect to VR(X) yields the projection operator 

E ($ Q M Q ^ <5 2 )=E (E~ =1 (6/n 2 a 2 ) $ 2 < 5 2 ) = E ($ 2 < 5 2 ) , 

as desired. This is an (extreme) example of how the present procedure deals 
with a reducible set of constraints without needing to extract an irreducible 
subset. For more complicated examples it may be useful to experiment with 
different choices of M alS . 



The full-time physical evolution operator 

Up to this point we have determined how to construct the desired projection 
operator E for any (nonzero) short-time interval e = T/N. Let us assume 
that we have carried out the necessary Lagrange multiplier variable inte- 
grations for every time slice, i.e., all n, 1 < n < N. The result will be 
the evolution operator restricted to the physical Hilbert space — let us call it 
S(T) — which is given by 

S(T) = lim Te^-C-D^^E e - l C-Z n{t)dt E • • • E e ^^ n{t)dt E . 
The indicated limit can be taken with the result that 

£(T) = TEe^o (EM(1)E)(! 'E =E [T e~^l {Enm) dt ] E . 

This result holds whether or not the effective Hamiltonian He (t) = E7i(t)E 
at time t is (essentially) self adjoint or merely symmetric (Hermitian) on 
fjphys- If it is not self adjoint, the resulting evolution will not be unitary. On 
the other hand, there are many cases of interest when E7i(t)E will be self 
adjoint for alH, < t < T, and the evolution on the physical Hilbert space 
will be by unitary operators £(T) acting on f) p h ys , and which, in fact, act 
as the zero operator on the Hilbert space S) p hys = i]0 -ftphys, namely, the 
subspace othogonal to fjphys, which may be called the "unphysical" Hilbert 
space i^unphys = -ftphys- This follows from the fact that Efjp hys = 0. For con- 
venience, we will assume that the effective Hamiltonian is self adjoint and 
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therefore that the resultant evolution is unitary. (Note, even when the effec- 
tive Hamiltonian is not self adjoint, that fact alone does not imply that the 
resultant evolution takes physical vectors into unphysical vectors; whatever 
the form of evolution, vectors that start in the physical Hilbert space will 
remain in the physical Hilbert space.) 



Constructive approach 

It is very instructive to see the lattice-limit construction for the evolution 
operator combined with the constraint operators in as direct a form as pos- 
sible. For this purpose we assume that Ti.(t) is a continuous function of t. In 
that case, and modulo domain considerations, we assert that 

E [Te-Cf™^*]! 

= l }™ [ II dJ{r n )dT{ ln ) lim I J] dS-(A r(n , ) 

6_> l<n<N 6 ~* J 1<1<N' 

x[l - ie"H r{n:l) - ie'\ a r{n ^ a - -l^K^Xin^P +■■■] 

Here the single index n on the variables r and 7 labels consecutive time steps 
with a spacing e. Observe that only the first-order term (i.e., proportional 
to e') is necessary to include for the Hamiltonian, but to properly include 
the constraints it is necessary to include first- and second-order terms (i.e., 
proportional to e' and e' 2 ). This difference between the dynamics and the 
constraints holds because the integration over the Lagrange multiplier vari- 
ables is, in reality, an integral over a (regulated) ^-correlated distribution. 

This formula may also be brought to yet another interesting form. By 
carrying out the integrations over all \ r ( n ,i) and 7„ we are left with the ex- 
pression 

E [Te"^ (EH(t)E)d *]E 

= lim f Y[ e-iMn-ieM^M-^p) dJ ^ ^ 

J l<n<N 

where H n = ?i r r n ,i) an d 



/< 



.)<n>„)= lim lim f\...) ^l±(hi iTn . 

C^0+ J^ L 7TT n 
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This expression may be a good place to start calculations in certain circum- 
stances. 

Special case 

A great deal of simplification takes place whenever 7i(t)K = K7i(t) for all t, 
< t < T. In this case we have the identity 

E [Te-Ct^^^jE = [Te^ w W dt ]E . 

Observe that the evolution is then automatically unitary in the physical 
Hilbert space, and that it is only necessary to impose the projection operator 
E once, say, as we have done, at the beginning of the evolution, very much in 
the manner of an initial- value equation used in the case of classical first-class 
constraints. 

Remark 

We hasten to add that for a given, specific set of constraints it is generally 
possible to construct E by alternative and possibly simpler integral represen- 
tations. By all means, when simpler procedures exist they should be used. 
For example, as noted earlier in the case of a set of constraint operators which 
close to form a Lie algebra for a compact group, a group invariant average 
over all unitary transformations continuously connected to the identity and 
generated by the constraint operators leads to the same projection operator 
E given by the general prescription developed in this paper for sufficiently 
small 5, 5 > 0. 

Relation with Dirac's prescription 

The machinery developed above is expressed entirely in an abstract operator 
form plus a functional integral over the Lagrange multiplier variables which 
has a well-defined formulation as a lattice limit. Another abstract approach, 
and one of the more common methods for dealing with constraints in the 
presence of dynamics, is that proposed by Dirac ||. Briefly summarized, 
Dirac's procedure defines the physical Hilbert space as the subspace spanned 
by vectors \ip) p h ys that fulfill the requirement 

|^)phys = 
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for all a, implying an exact fulfillment of the quantum constraints. For 
second-class constraints, such as $1 = P and $2 = Q, where P and Q 
are conventional Heisenberg variables, there is no nontrivial solution of the 
equations P|?) = Q\?) = for that would require [Q,P}\7) = ih\?) = 0. 
Moreover, for a constraint $ with zero in its continuous spectrum, the equa- 
tion $|#) =0 has no (normalizable) solution in the Hilbert space, i.e., there 
is no nontrivial solution for \ jf). Thus, for second-class constraints, and for 
general, open first-class constraints, as well as for constraints with zero in the 
continuum, the basic Dirac prescription is unsatisfactory. Of course, there 
are modifications of the Dirac scheme, such as the use of a Dirac bracket 
[||] in the classical theory followed by the adoption of this bracket in the 
process of canonical quantization. Unfortunately, the classical phase space 
associated with the Dirac bracket generally has a nontrivial curvature (and 
possibly a nontrivial topology), and so a straightforward application of the 
quantization rules generally leads to an incorrect answer. 

Nevertheless, there are some elements of the Dirac procedure and the 
procedure developed in the present paper that are rather close to each other, 
at least in spirit. For one thing, we note that if $ Q |?) = for all a, then 
($Q,£ a/3 $/3)|?) = and conversely, where summation is implied and we are 
only considering the case M al3 = 5 af3 for simplicity. For the second-class 
constraints considered above we would be asking that (P 2 + Q 2 )|?) = 0, 
which, as is well known, implies that |?) = 0, an unsatisfactory solution. 
Let us relax this condition (cf., f|) and say that (P 2 + Q 2 )|?) = for 
which we do have a nonvanishing solution, namely |?) = |0), the ground 
state of the harmonic oscillator. There is no violation of the classical limit 
for these constraints by this modification since, as h — > 0, we recover the 
original, classical, second-class constraints, namely, p = and q = 0, in a 
standard notation. However, the generalization of the Dirac procedure just 
illustrated suffers from one important defect: The extended criterion does 
not satisfactorily specify the physical Hilbert space. Take the given example 
and instead let (P 2 + Q 2 )|?) = ch\l) for some constant c. If c = 3 there 
is a solution which is |?) = |1), the first excited state, but for c = 2 — or 
indeed for any c, 1 < c < 3 — there is no nontrivial solution and so |?) = 0, 
even though for any of these cases the recovery of the classical second-class 
constraints still applies. 

Rather than relax the original Dirac prescription in the manner just dis- 
cused, we have chosen to look for solutions to the equation E($ Q 5 a/3 $/3 < 
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5 2 (fo))|?) = I?). It is convenient to examine this expression as a function of 5. 
For sufficiently large S, there is generally a nontrivial subspace of vectors sat- 
isfying this relation. As 5 decreases, the space in question remains the same 
or its span may decrease. As an example, again take the second-class case 
and consider E (P 2 + Q 2 < ch)\l) — |?). Beginning with small c, we observe 
for c < 1 that there is no nontrivial solution, |?) =0; for 1 < c < 3, there is 
a single nontrivial solution, |?) = |0), which spans a one-dimensional Hilbert 
space; while for 3 < c < 5, there is a two-dimensional Hilbert space spanned 
by the two vectors |0) and |1), where the latter vector is the first excited 
state of the harmonic oscillator; etc. This dependence on the parameter c is 
to be preferred to that described above. 

Suppose the second-class constraints are irregular as well. In particular, 
let $! = P and $ 2 = Q 3 . We are then led to consider E (P 2 + Q 6 < 
c'tfl 2 ). Let c' = c' denote the value that captures the ground state of 
the "Hamiltonian" P 2 + Q 6 , let d = c[ denote the value that captures the 
first excited state, etc. The argument then proceeds as above. Of course, 
E(P 2 + Q 6 < c'h 3 / 2 ) ^ E(P 2 + Q 2 < ch), for general c and c', but this 
fact merely represents the ever-present 0(h) ambiguity in any quantization 
procedure. 

Time-dependent constraints 

In the time-dependent case $ Q — > for all a and t, < t < T, and 

assuming each constraint operator is continuous in t as well as the needed 
self-adjoint properties, we are quickly led to an expression for the evolution 
operator in the time- dependent physical Hilbert space, given by 

£'{T) 

= lim Te" 1 j C"-D« W(t) dt E N -i e'^ >« W(t) dt E N _ 2 ■ ■ ■ E , e~ l J ^ dt E . 
Here 

E n EEE($ Q (ne)M^(ne) < 6 2 ) 

/r (n + l)e 
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While this expression is correct, there does not seem to be any simpler and 
more elegant operator expression that holds in the general case of time- 
dependent constraints. 

We do not pursue time-dependent constraints further, and hereafter we 
again assume that all the constraints have no explicit time dependence. 

3 Functional Integral Representation 

It is convenient to adopt a functional representation for the original Hilbert 
space $), and to present the main properties of the foregoing argument in 
that language. For this purpose, let us focus on a canonical system with 
dim(i}) = oo and J degrees of freedom, J < oo, described by the basic, 
irreducible, self-adjoint operators P = {P J }and Q = {Q-*}, 1 < j < J, which 
satisfy the canonical commutation relations [Qi , P k ] = i8^ k \ with all other 
commutation relations vanishing. As an irreducible set of operators, all other 
operators may be constructed from the basic ones. Hence, TC(t) = TC(P, Q, t) 
and $ a = Q) for some suitable functions. 

Although other bases may be used, we concentrate here on the functional 
representation induced by coherent states. Thus, along with the canonical 
operators, we introduce canonical coherent states defined by 

\p,q)=e-^ p e^ Q \ V ), 

where p = {p 7 }^ =1 and q = {q^}j =1 , with each p 7 and g J an arbitrary real 
number, q-P = S^ =1 q*P\ etc., and \r]) is a fixed, normalized fiducial vector. 
Moreover, for any choice of \rj) it follows || that 

/MNP, 9 |«W*,) = i, W»,) = *U#*tl2*. 

Needed domain restrictions on \if) will be assumed, without comment, as 
they arise. At this point p and q, although having suggestive names, are 
just mathematical labels of the coherent states. As a first step in giving 
meaning to p and q, we assume that \rj) is "physically centered", meaning 
that (?7|P|?7) = and (r)\Q\r)) = 0. In this case it follows that (p, q\P\p,q) = p 
and (p,q\Q\p,q) — q- As a second step, we assume that \rj) is "physically 
scaled", i.e., that (r]\P 2 \ri) = 0{h a ) and (v\Q 2 \v) = 0{h h ), where a > 0, 
b > 0, and a + b = 2. 
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To obtain the desired functional representation for the lattice formulation 
we repeatedly insert coherent-state resolutions of unity leading to 



In this expression we have again used r(n,l) = (n — 1)N' + (I — 1) and 
identified (p",q") = (pnn',Qnn') an d (p',q') = {Po,qo}- We reemphasize that 
it is necessary to expand the constraint operators to second order in the 
lattice spacing in order to correctly obtain the desired projection operator. 

It is customary to interchange the order of integration and the continuum 
limit in order to arrive at a formal expression for the propagator in the 
presence of the constraints. In so doing one adopts for the integrand the 
expression it would assume for continuous and differential paths (even if 
such paths do not contribute in the continuum limit!). In so doing one is 
led to identify the classical action for which such paths are appropriate. If 
we proceed in the same fashion, we arrive at the formal functional integral 
expression 



We further note that H{p,q,t) = H(p,q,t) + 0{ff), c > 0, and <fi a (p,q) = 
®a{P) q) + 0(h d ), d > 0, thus leading to the expected behavior as h — > 0. 

It is clear from this expression that the system we have quantized here is 
determined by the (H augmented) classical action functional 



(p",q"\E[Te' t fl (Enm)dt ]E\p',q'} 





In this expression we have introduced the symbols 



H(p,q,t) = (p,q\H(P,Q,t)\p,q) 

<t> a {p,q) = {p,q\$a(P,Q)\p,q} ■ 



Jo \P-Q - H (P> 9, t) - X a (t) <f) a (p, q)] dt , 
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which has the expected form of a canonical system in the presence of a 
number of constraints. Evoking a stationary action principle for this classical 
action leads to both dynamical and constraint equations as usual. Note 
that the second-order terms in the constraint operators needed in the lattice 
formulation to construct the correct quantum projection operator disappear 
in this formal expression since, for paths and Lagrange multiplier variables 
that are smooth and differentiable, i.e., have classical behavior, these second- 
order terms formally disappear. 

While the derived classical action essentially reflects the original quantum 
theory, it is necessary, as always, to exercise great care in elevating a formal 
path integral involving the classical action to a valid calculational procedure. 
While not every lattice-limit formulation will lead to the correct answer, we 
emphasize that one acceptable construction procedure is the one presented 
in this paper. 

Commentary 

In a classical theory, such as the one above, it is true that the Lagrange multi- 
plier variables play a significant role in the solution of the classical equations 
of motion, either being fixed functions chosen to enforce the constraints in the 
case of second-class constraints, or needing some arbitrary specification (se- 
lection of a gauge) when the Lagrange multiplier variables are not determined 
by the equations of motion in the case of first-class constraints. However, 
what holds true in the classical theory has no reason to hold true in the 
completely different quantum theory. For example, it should be noted that 
nowhere in our discussion have we used the words "gauge" or "gauge fixing" 
in achieving the proper quantum results. Indeed, in the abstract formulation 
of Section 2, the concepts of "gauge" and "gauge fixing" are completely for- 
eign to the whole procedure. If anything, "gauge fixing", and its associated 
difficulties, is a classical concept inappropriately and unnecessarily carried 
over into the quantum theory. In concrete terms, and with respect to the 
formal phase space path integral we have derived, we note that no discussion 
of gauge invariance of the formal path integral has been given nor is any 
needed. 

If the constraints are all first class and closed, then gauge invariance of 
the resulting theory as expressed in the (true) physical Hilbert space may 
be demonstrated |]. On the one hand, if 5 > is chosen sufficiently small 
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so that E ($ a M a P$p < 5 2 ) = E ($ a M al3 $p = 0) ^ 0, then invariance of 
the resulting theory in the physical Hilbert space follows from the fact that 
e -itv a <s> a £ _ £ £ Qr arbitrary real-valued {c<J a }^ =1 . On the other hand, if 
E($ Q M Q/3 $ /3 < <5 2 ) ^ Ef^^^ = 0) no matter how small 5 > is, 
then it means that e _lwa * a E = E + 0(5 -E), and gauge invariance arises 
in the true physical Hilbert space i)trucph y s only after a suitable 5^0 
limit as illustrated in Appendix B. Of course, for second-class constraints 
the issue of gauge invariance does not arise. The presence of such constraints 
is distinguished by the fact that E = for all 5 < 5' for some 5' > 0. 

Additionally, we observe that the abstract operator formalism developed 
in Section 2 applies to a field theory just as well as to general systems of 
finitely many degrees of freedom, canonical or otherwise. To discuss any 
system — a model field theory included — one only needs to have candidate 
operators for the Hamiltonian and for the various constraints expressed in 
the original Hilbert space. 
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Appendix A 



The construction of the first projection operator example in Section 1 relies 
on a conditionally convergent integral, defined for — oo < x < oo, which is 
given by 

We may assign one meaning to this expression by choosing the formula 

hm f L e-^dA, 

an expression which defines the function 

fi(x) = 1 if \x\ < 5 , 
= | if \x\ — 5 , 
= if \x\ >5 . 

It is clear that fi(x) 2 ^ /i(^), specifically at \x\ = 5, which would be needed 
to ensure that E 2 = E . 

To overcome this difficulty, we adopt the prescription 

r / \ v v f L -ix\ sin[(5 + C)A] 
fi(x) = hm hm / e lxA — - — , - dX 

C^0+ L^oo J_ L 7rA 

= 1 if \x\ < 6 , 
= if \x\ > 5 . 

Clearly in this case f\(x) 2 = f\(x) for all x. Likewise we can also define 

fo(x) = 1 if \x\ < 5 , 
= if |x| > £ , 

which also satisfies /o(^) 2 = fo( x ) for all x, by replacing lim^^ + by lim^o- 
in the definition of f\{x). 

In application to a single, self-adjoint quantum constraint $, we conclude, 
for any 5 > 0, that 

E($ 2 <5 2 )= lim hm / L e-^ Sin[( ^| C)A] g?A 

C^0+ L^oo J_ L 7rA 
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and 

E ($ 2 < 5 2 ) = lim lim I" C)A] dX . 

In many cases the resulting projection operators are identical. An example 
for which they are not identical is given by 

|0) (0| = E (P 2 + Q 2 < ti) j^E (P 2 + Q 2 < h) = . 

Here the one-dimensional projection operator |0)(0| projects onto the ground 
state for a harmonic oscillator "Hamiltonian" P 2 + Q 2 . 

Unless stated otherwise, we shall adopt the convention that 

r e ^^dA,E($ 2 <5 2 ). 

J-oo n * 

Appendix B 

The issue of dealing with constraint operators with their zero in the contin- 
uous spectrum has been dealt with by several authors |§. Here we offer a 
brief discussion in keeping with the style of the present paper. 

Let us focus on the projection operator E (X 2 < 5 2 ) where X = lies in 
the continuous spectrum of the self-adjoint operator X defined in the original 
Hilbert space f). We write 5j = fh <8> S) 2 such that, for sufficiently small S, 
X = Xi (8> l2> where Xi is nondegenerate while S)2 accounts for possible 
degeneracy. Generally, this factorization depends on 5, but we assume that 
a suitable limit of these spaces exists as 5 — > 0. In the same decomposition, 
E = Ei ® 1 2 . For product vectors = \ipi) <8> \ip2) £ it follows that 

(^|E(X 2 < 5 2 )\<P) = (^ilE^X 2 < 5 2 )l0iX^|02> 

J-S 

where E\(x) is the spectral family for X%. For the first factor we choose a 
functional representation that diagonalizes X\ leading to 




ipi(x)*<pi(x) dx . 
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For every j^i) and \<f>i), this expression vanishes as 5 — > 0. However, before 
taking the limit let us first rescale this integral so that it reads 




ipi(x)* (fii(x) dx . 



This expression has a well-defined and generally nonvanishing limit for a sub- 
set of elements ipi(x) and (f>i(x), namely, those that are continuous functions 
in an open set about x — 0. For the subset of continuous functions, it follows 
that ^ 

limi- / MxTMx)dx = MoyMo) ■ 

s^o 2d J_s 

To summarize, we determine, for suitable elements j^i) and \(j>i), that 
lim(25)- 1 (^|E(X 2 < 5 2 )|0) = Vi(0)>i(0) (V> 2 |02> • 

If this expression is regarded as an incomplete inner product, then, after 
completion, we learn that the resultant space is C ® i} 2 — #2 since C, the 
space of complex numbers, has only one (complex) dimension. We conclude 

that f)truephys = -£)2- 

We can understand this result in more physical terms. Let us realize 
Sji as L 2 {\— 5, 5]) and introduce an orthonormal basis given by u n (x) = 
(1/V25) e inwx/s , n e {0, ±1, ±2, . . .}. If the Hamiltonian contained the term 
P 2 = —d 2 /dx 2 , or such a term was added to the Hamiltonian, then 

P 2 u n (x) = (nix / 5) 2 u n (x) 

since u n is an eigenfunction of P 2 . For very small 5, all energy contributions 
are huge save for the single case n = when this contribution vanishes. 
Thus, when <5 <^ 1 energetics favor only the ground state u (x) = l/\^26 in 
an expansion of ip\(x) and (f>i(x), and consequently 

= ^i(O)*0i(O) 

for functions which are continuous about x — 0. 
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Is it really necessary to take the limit 5 — » 0? In our view the practical 
answer is "No". For example, if energies are restricted to those which are 
much less than 5~ 2 , a number which can be chosen as large as one likes, then 
only the ground-state contribution arises for S)i and hence i^truephys = -fjph yS! 
as desired. 
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